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Abstract 


A  Monte  Carlo  Fokker  Planck  Code  is  used  to  investigate  the 

scattering  loss  of  particles  from  a  relativistic  electron  beam 

propagating  in  a  low-pressure  air  channel.  Two  special  cases  are 

examined.  The  first  corresponds  to  an  experiment  in  which  a  one- 

cm  radius,  90  keV  beam  propagates  in  a  hydrogen  channel  of  the 

14  -3 

same  radius,  50  cm  long  with  a  density  of  10  cm  and  a 
temperature  of  0.5  eV  assuming  the  interaction  to  take  place  only 
with  the  charged  particles  of  the  medium.  It  is  shown  that  the 
beam  will  not  traverse  such  a  channel  if  the  ionization  drops 
below  20%  noting  that  for  such  parameters  the  channel  pressure 
corresponds  to  few  tens  of  militors.  The  second  case  corresponds 
to  what  is  referred  to  as  the  Ion  Focused  Regime  (IFR)  where  the 
pressure  is  about  0.1  torr.  Scattering  of  beam  particles  by  both 
neutral  and  charged  particles  is  taken  into  account  in  this  case 
as  well  as  multiple  scattering  from  neutral  targets.  It  is  shown 
that  transmitted  fraction  is  linearly  proportional  to  the  square 
of  the  channel  radius,  and  the  probability  of  multiple  scattering 
increases  with  increasing  channel  radius. 

Introduction  and  Analys i s 

In  many  applications  involving  the  use  of  charged  particle 
beams  propagating  in  a  gaseous  medium  the  question  often  arises 
as  to  the  efficiency  of  propagation  and  the  processes  that 
detract  from  efficient  transmission.  Of  special  interest  in  this 
connection  is  the  propagation  of  intense  relativistic  electron 
beams  in  low-pressure,  pre-ionized  media.  Although  there  are 


several  mechanisms  that  impact  such  a  propagation  the  effect  of 
scattering  is  especially  important  and  requires  extensive 
investigation.  This  effort  is  aimed  at  studying  the  role  of 
scattering  in  electron  beam  propagation  in  low  pressure  channels. 
It  is  carried  out  with  the  aid  of  a  Monte  Carlo  Fokker  Planck^ 
(MCFP)  computer  Code  which  allows  us  to  follow  a  given  sample  of 
non-interacting  particles  in  full  configuration  and  velocity 
space  as  the  beam  tranverses  a  cyl indrical ly  symmetric  gaseous 
channel.  The  four  main  parts  of  (MCFP)  are:  The  Fokker  Planck 
section,  the  large-angle  scattering  section,  the  neutral 
scattering  section,  and  the  particle  loss  section. 

In  the  first  section  the  code  follows  a  sample  of  particles 
in  time  by  advancing  the  position  and  velocity  in  accordance  with 
the  collisionless  equations  of  motion.  Particle  deflections  and 
energy  depostion  in  the  background  plasma  are  statistically 
determined  using  the  usual  scale  times  derived  from  Fokker-Planck 
theory.  For  the  deflection  of  beam  particles  MCFP  assumes  that 
each  small  angle  scattering  results  in  a  deflection  angle  9 
(usually  set  to  5°  as  input)  and  determines  the  number  of 
deflections  encountered  by  the  equation: 
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where  &  t  is  a  given  time  step, 
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familiar  deflection  given  in  MKS  units  by 
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In  the  above  equation  tylx)  denotes  the  error  function  and 
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m  =  mass  of  test  particle 

m^  =  mass  of  particle  of  species  "b" 

=  temperature  of  particles  "b" 

Ze  =  charge  of  test  particle 
E  =  energy  of  test  particle 

/\D  =  Debye  length  = 

<rQ  =  permittivity  of  free  space 

The  above  small  angle  scattering  information  is  then  used  to 
update  the  velocity  vector  of  a  given  sample  particle.  The 
energy  loss  of  a  given  test  (beam)  particle  is  obtained  in  a  way 
similar  to  that  of  small  angle  scattering  just  described.  In 
this  case  if  we  let  E  be  the  energy  lost  by  a  sample  particle  in 
time  &  t  then 
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where  E  is  the  beginning  energy,  and  is  the  familiar  energy 


exchange  time  given  by 
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with  the  parameters  defined  as  before.  Thus  the  speed  of  the 
test  particle  can  also  be  updated  in  each  time  step.  It  might  be 
noted  at  this  point  that  MCFP  has  an  option  that  determines  the 
background  plasma  temperature  by  employing  the  simple  model  in 
which  the  channel  temperature  difference  (T-T  )  is  directly 
proportional  to  the  total  energy  deposited,  AEfc,  at  time  t. 

This  clearly  implies  that  the  background  plasma  instantly 
thermalizes  once  energy  is  deposited.  The  newly  determined 
temperature,  T,  is  then  used  to  calculate  the  new  ^  and  q 
assuming  that  the  background  plasma  remains  Maxwellian  and  thus  a 
beam-channel  feedback  effect  can  be  studied.  If  the  changing 
background  temperature  model  is  not  selected  then  MCFP  assumes  a 
constant  temperature  for  an  unperturbed  Maxwellian  background. 

The  code  also  takes  info  account  the  possibility  that  a 
beam  particle  suffers  a  l=>rge  angle  scattering  with  the 
background  plasma.  Usirm  a  random  number  the  code  determines 
such  an  event  by  samplin:  from  the  distribution. 
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event  in  the  center  of  mass  frame  is  then  obtained  by  randomly 


sampling  the  equation 


27T 


where 

<VeJ 


J  0"  f  kv  9' &  C*P 


(6) 

(7) 


and  is  the  Rutherford  cross  section.  The  quantity  0- 


represents  the  angle  that  defines  the  onset  of  large  angle 
scattering,  and  E  is  the  relative  kinetic  energy.  In  the 
interest  of  computation  any  deflection  angle  in  the  laboratory 
frame  less  than  10°  is  ignored.  When  the  kinetic  energies  of  the 
interacting  particles  is  above  a  certain  threshold,  Efch,  the  code 
treats  the  background  particle  as  a  test  (beam)  particle  meaning 
that  its  spatial  and  velocity  coordinates  are  tracked  in  time. 

If  the  above  requirement  on  the  kinetic  energies  is  not  met  then 
MCFP  will  relegate  the  energy  of  the  less  energetic  particle  to 
the  background  plasma.  After  that,  only  the  higher  energy 
particle  from  the  interaction  is  followed  allowing  us  a  method  by 
which  "delta"  rays  can  be  identified  and  followed. 

In  the  case  of  scattering  from  neutral  particles  the 
probability  for  such  an  interaction  is  selected  by  first  defining 
an  important  mean  free  time,  namely 
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where  \J  is  the  beam  electron  velocity,  n  is  the  neutral  density 
andT is  an  appropriate  interaction  cross  section.  Although  other 


inelastic  events  could  be  considered  we  limit  our  study  to 
ionization  and  the  interaction  probability  can  then  be  described 


by  dt/tf  when  dt  is  taken  as  a  very  small  code  time  step.  Given 
an  appropriate  differential  scattering  cross  section  a  random 
scattering  angle  ©  can  be  selected  from  the  distribution 
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where  1[  is  a  random  number.  For  relativistic  electrons 
propagating  in  a  low  temperature  medium  it  is  expected  that  the 
distribution  be  highly  peaked  in  the  forward  direction,  and  using 


the  results  from  hydrogen  appropriately  sealed  for  higher 
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elements  we  can  write 
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where  in  CGS  units  the  parameters  £  and  (3  are  given  by 
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In  the  above  expression  m  is  the  electron  mass,  k  is  the 
Debroglie  wave  number,  h  is  Planck's  constant,  Q~  —  /•  2 

is  the  Bohr  radium,  and  is  the  minimum  scattering  angle 

given  by 
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where Y  is  the  relativistic  parameter  and  c  is  the  speed  of 
light.  As  we  noted  earlier,  small  angle  scattering  is  included 
in  the  code  by  an  implicit  Fokker-Planck  method  where  the 
probability  for  such  a  scattering  event  in  the  code  time  t  is 
simply  where  td  is  given  by  Eq.  (2),  for  plasma  densities  on 

M 

the  order  of  1011  -  1013  cm  3  which  might  be  appropriate  for  IFR 
and  for  l/~C  it  is  readily  seen  thatTd^lO  ^  seconds  so  that  at 
pulse  lengths  on  the  order  of  nanoseconds  it  is  clear  that 
scattering  from  charged  particles  can  be  neglected.  These 
considerations  are  consistent  with  some  experimental 
observations .  ^ 3  ^ 

In  order  to  crudely  model  the  propagation  of  a  relativistic 

{ 4 ) 

electron  beam  in  either  an  IFR  or  a  force-free  propagation 
mode  we  have  incorporated  the  self  magnetic  field  in  a 
phenomonological  way  using  Ampere's  law  for  an  infinite,  uniform 
beam.  We  have  neglected  the  radial  electric  field  by  taking  the 
charge  neutralization  factor  f  to  be  unity.  Since  the  code 
calculates  the  number  of  particles  that  traverse  the  channel,  the 
fraction  of  the  beam  that  is  transmitted  would  then  represent  an 
upper  limit  for  pulse  lengths  that  are  shorter  than  the 

charge  neutralization  time,  and  a  good  estimate  when  since 

the  instabilities  which  might  occur  in  such  time  scale  are 
ignored.  Under  these  conditions  the  description  of  the 
scattering  process  might  be  deemed  appropriate  to  that  taking 
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place  in  a  force  neutralized,  charge  neutralized  or  force-free 
regimes.  With  propagation  having  been  observed  in  pressure 
ranges  of  0.1  -  5  torr ^  it  should  be  noted  that  at  lower 
pressures  there  are  usually  insufficient  neutrals  to  provide  an 
adequate  ion  density  for  force  neutralization  and  at  higher 
pressures  multiple  scattering  could  cause  dispersal  of  the  beam. 
The  latter  effect  has  been  incorporated  in  the  code. 

Finally,  it  should  be  noted  that  in  assessing  the 
usefulness  of  MCFP  code  in  beam  propagation  studies  it  should  be 
kept  in  mind  that  particle  scattering  can  be  accurately  described 
and  that  no  assumptions  need  be  made  concerning  the  plasma 
conductivity  since  ionization  can  be  calculated  self- 
consistently .  Particle  and  energy  loss  mechanisms  can  also  be 
readily  included  so  long  as  the  appropriate  mean  free  paths  for 
these  processes  can  be  selected.  Because  of  its  three 
dimensional  simulation  capability,  the  often  used  "paraxial” 
approximation  and  others  need  not  be  made.  Because  of 
limitations  on  the  number  of  particles  that  can  be  simulated  it 
is  clear  that  the  treatment  of  self  fields  and  the  modelling  of 
instabilities,  though  possible,  might  be  quite  difficult. 

Some  Appl icat ions 

When  used  in  connection  with  the  study  of  an  electron  beam 
propagation  in  a  low  pressure,  preformed  channel  given  by  the 
following  parameters: 

beam  current  =  lkA 

beam  radius  =  1  cm 
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beam  energy  =  90  keV 
plasma  species  =  hydrogen 
channel  radius  =  1  cm 
channel  length  =  50  cm 

14  -2 

channel  proton  density  =  10  cm 

14  -3 

channel  neutral  density  =  4  x  10  cm 
channel  temperature  =  0.5eV 

__  g 

pulse  length  =  0.8  x  10  sec 
The  code  revealed  that  no  propagation  occurs  below  20% 
ionization.  At  the  above  parameters  the  results  show  that  about 
33%  of  the  beam  particles  were  last  in  a  distance  of  about  10cm 
down  stream  from  the  injection  point.  When  the  neutral  density 
was  reduced  by  a  factor  of  10  most  of  the  beam  particles  were 
able  to  traverse  the  channel. 

In  assessing  the  rate  of  multiple  scattering  by  neutral 
particles  the  following  case  was  examined: 
channel  radius  = 

beam  radius  H  =  1.5  cm 

b 

neutral  density  Y1  =  2  x  lO^cm  ^ 

n 

electrons  used  to  simulate  the  beam  H  =  5000 

s 

Transmitted  fraction  of  beam  =  T 

Number  of  scatter::  as  off  of  a  neutral  particle  = 

The  ratio 

The  results  are  a;  ;r.  Table  1. 
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Table  1 


K 

T 

m 

S 

2.5 

0.007 

7338 

1.47 

5 

0.015 

9258 

1.85 

10 

0.03 

12183 

2.4 

15 

0.079 

15003 

3.0 

21.2 

0.121 

17807 

3 .56 

We  note  that  the  transmitted  fraction  is  almost  linearly 

proportional  to  which  is  a  direct  consequence  of  the  geometric 

effect.  Another  geometric  consequence  is  the  linear  scaling  of 

and  S  with  A.c  .  Since  the  probability  of  neutral -interaction  is 

proportional  to  the  path  length  traversed,  and  this  length 

increases  with it  is  expected  that  multiple  scattering  should 

increase  with  .  The  variation  of  these  parameters  with  the 

neutral  density  is  shown  in  Table  2  for  A.  =  21.2  and  A.,  =  1.5 

1  c  b 

cm  and  =  3  x  10  ^cra^ 

tot 

Table  2 


(cm-3 ) 

l 

m 

S 

T 

2x1013 

293 

KO 

o 

o 

0.95 

2x1014 

2789 

0.66 

0.79 

2x1015 

18000 

3 . 6 

0.14 

2x1016 

68263 

13.65 

0.002 
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